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ABSTRACT 


Like  any  other  complex  dynamic  system  the  human  body  responds  in  a  complex 
way  to  acceleration  inputs  which  vary  rapidly  with  time.  The  need  to  avoid  stresses 
large  enough  to  cause  injury  to  the  body  usually  imposes  limits  on  the  permissible 
input  acceleration. 

The  restraint  system  interposed  between  a  vehicle  and  its  occupant  can  modify  the 
physiological  effects  of  a  vehicle's  acceleration  -  time  history.  This  modification 
should  be  made  as  favorable  as  possible  by  minimizing  the  etresses  generated  in 
the  vehicle's  occupant.  To  determine  optimum  dynamic  characteristics  for  the  re¬ 
straint  system,  its  important  characteristics,  and  those  of  the  human  body,  need  to 
be  represented  in  terms  of  a  mathematical  or  "d3mamic"  model.  Through  suitable 
analysis,  either  mathematical  or  by  means  of  a  computer,  those  dynamic  charac¬ 
teristics  of  the  restraint  system  can  be  determined  which  will  minimize  the  peak 
stresses  developed  in  its  human  occupant. 

hi  this  report  a  general  theory  of  suitable  dynamic  models  is  developed  for  th's 
type  of  problem.  Closed  form  solutions  for  a  number  of  simple  cases  are  pre¬ 
sented  also.  In  addition  a  method  is  shown  which  permits  development  of  simple 
dynamic  models  for  the  human  body  utilizing  existing  e3q)erimental  data. 

Most  test  data  has  limitations.  This  seems  particularly  true  when  the  subject  is 
as  complex  and  vai^iable  as  the  human  body.  The  limitations  associated  with  the 
application  of  physiological  data  to  dynamic  models  of  the  human  body  can  be  min¬ 
imized, however,  if  the  test  program  is  designed  with  this  application  in  mind. 
Accordingly,  the  reader  will  find  discussed  in  some  detail,  the  necessary  re¬ 
quirements  for  short  period  acceleration  testing  with  live  human  subjects  as  well 
as  some  suggested  requirements  for  dynamic  characteristics  of  test  rigs. 
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SECTION  I 


INTRODUCTION 


This  investigation  is  concerned  generally  with  the  stresses  developed  in  the  body 
of  a  human  occupant  of  an  aerospace  vehicle  which  is  subjected  to  short-period 
acceleration.  More  specifically,  its  objective  is  to  ensure  that  the  characteristics 
of  the  ’’restraint  system,"  interposed  between  the  occupant’s  body  and  the  vehicle, 
will  minimize  body  stresses. 

To  accomplish  this,  a  means  has  been  devised  of  calculating  the  response  of  the 
body  and  associated  restraint  system  to  the  applied  acceleration.  A  mathematical 
model  has  been  developed  which  represents  those  characteristics  of  the  body  and 
restraint  system  which  are  significant  in  the  study  of  the  problem.  Such  a  model 
necessarily  involves  simplification  and  represents  only  the  major  system  consid¬ 
erations.  Nevertheless,  it  has  considerable  value  in  predicting  limiting  stresses; 
indeed,  the  variability  of  human  bodies  and  restraint  system  materials  >  may  not 
justify  the  use  of  a  more  complex  model. 

The  most  elementary  dynamic  system  in  the  field  of  restrained  body  dynamics, 
illustrated  in  figure  l,is  mathematically  quite  sophisticated.  Although  partially 
closed  form  solutions  can  be  obtained  for  a  great  many  cases,  their  derivation  is 
often  so  lengthy,  and  the  resulting  equations  so  long,  that  the  effort  involved  is 
hardly  worthwhile  in  an  applied  research  program  where  rapid  results  are  required. 


SYSTEM  (2) 


SYSTEM  (1) 


Figure  1.  Elementary  dynamic  model  in  series 
with  a  restraint  system. 


For  this  reason,  the  influence  of  a  particular  restraint  system  is  usually  investi¬ 
gated  with  the  aid  of  a  digital  or  analog  computer.  Despite  these  aids  of  modern 
technology,  it  is  still  necessary  to  organize  the  basic  equations  of  motion  if  we  are 
to  devise  digital  or  analog  logic  of  maximum  economy.  It  is  also  necessary  to  de- 


rive  closed  foiiu  solutions  for  limit  cuses  ui  order  tust  the  results  of  computer 
runs  can  be  checked  during  the  "de-bugging"  process,  and  to  provide  general  de¬ 
sign  "guide-lines"  as  to  the  most  promising  type  of  system. 


The  basic  djmamic  equations  involved  in  the  problem  were  identified.  From  these 
equations,  closed  form  solutions  for  certain  limit  cases  were  derived  to  permit 
checking  ot  computer  solutions.  Concurrently,  the  se  test  solutions  give  an  insight 
into  the  basic  physical  behavior  of  the  system. 


No  attempt  is  made  here  to  ejiplain  the  physical  implications  of  the  theory  to  the 
nonmathematician,  since  this  has  been  done  effectively  elsewhere. 


Simple  solutions  have  been  obtained  for  the  influence  of  slack  or  preloading  in  a 
restraint  S3rstem.  These  results  are  summarized  in  figures  11  and  13. 

By  definition,  a  restraint  system  degree  pf  freedom  consists  of  a  mass,  together 
with  a  single  spring  element  and  a  single  damping  element.  Both  of  these  rtiay  be 
discontinuous  and  nonlinear.  They  are  in  series  with  a  dynamic  model  which  rep¬ 
resents  the  human  body.  The  dynamic  model  may  have  any  number  of  degrees  of 
freedom,  both  in  parallel  and  in  series.  For  the  purposes  of  this  report,  only  the 
series  case  has  been  considered  in  detail  because  it  is  the  only  one  so  far  used  in 
practice. 

From  the  point  of  view  of  pure  mathematics,  there  are  two  types  of  restraint 
systems,  represented  by  continuous  and  discontinuous  functicas  respectively. 

When  the  function  is  continuous  and  differentiable,  one  set  of  differential  equations 
governs  the  entire  motion.  When  it  is  discontinuous,  two  different  sets  of  equations 
govern  the  motion  each  side  of  the  discontinuity,  and  in  addition,  the  change  from 
one  set  to  the  other  involves  generating  a  scries  of  initial  conditions  for  the  new 
equations. 


In  the  work  which  follows,  we  shall  first  investigate  the  dynamic  problem  in  a 
very  general  way,  imptosing  only  the  requirement  that  the  restraint  functions 
shall  be  real  and  single-valued.  Next  we  shall  consider  the  special  case  of 
linear  systems,  and  solve  a  number  of  simple  cases,  involving  both  linear  re¬ 
straint  with  discontinuities  (bottoming  systems)  and  continuous  nonlinear  restraint. 


While  it  is  fairly  easy  to  determine  the  dynamic  characteristics  of  a  restraint  sys¬ 
tem  by  physical  measurement,  the  human  body  presents  obvious  difficulties.  In 

X'axt  JLL  %JX  IrlUO  WC  XUOW'UOO  %fXlXO  ixuvr 

can  be  obtained  from  the  available  test  data.  Test  rig  criteria  for  future  e3q)eri- 
mental  programs  are  also  developed  from  theoretical  considerations ,  together  with 
the  theory  of  bounce  testing,  which  is  suggested  as  a  new  experimental  technique. 

It  is  felt  that  bounce  testing  offers  hope  of  measuring  bio-dynamic  characteristics 
which  are  not  easily  obtained  from  currently  employed  experimental  techniques. 
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SECTION  n 


DYNAMIC  THEORY 

1.  GENERAL  THEORY  OF  RESTRAINT  DYNAMICS 

a.  Generalized  equations  for  a  multidegree  of  freedom  system  in  series. 

In  body  dynamics  we  are  almost  always  concerned  with  dynamic  systems  in  series,  rather 
than  in  parallel.  The  springs  and  dampers  can  be  either  linear  or  nonlinear.  Such  a  sys¬ 
tem  has  a  number  of  characteristics  which  are  independent  of  the  characteristics  of  the 
springs  and  dampers.  We  shall  deal  with  these  before  proceeding  to  specific  relationships 
between  force,  deflection  and  velocity. 


For  the  system  in  figure  2,  let 

# 


value  of  the  mass  -A 
deflection  of  spring  </• 
velocity  of  spring  4^ 
force  in  spring 
force  in  damper 
5^'  the  maes  ratio. 


(3) 


(1)  ^ 


t 

d  ^ 

L...  1 

Figure  2. 


Series  dynamic  system. 


3^ 


3 


The  functions  -/y  and  ^  may  be  any  continuous,  single-valued  real  functions,  and  there 
need  be  no  relationship  between  the  functions  for  different  springs  and  dampers.  In  other 
words,  we  are  dealing  with  a  generic  system. 

If  we  now  write  the  equations  of  motion  for  each  system,  the  following  set  of  coupled  equa¬ 
tions  is  obtained: 

y,  -  -fA)  *  ^ 

Also 

Sr  ®  V  "■  ^  (2) 


where  is  the  unstretched  spring  length.  Differentiating  equation  2  twice  with  respect 
to  time, 


therefore 


#  # 


and  finally 


Substituting  for  in  equation  1 , 


4‘~\ 


I 


(3) 


(4) 


(5) 


or 


(6) 


•/  INPUT 


*  i 


ETC. 


Fif^ure  3.  Flow  chart  for  multidcpreu  of 
froedor>\  system. 


r: 


On  the  left  of  equation  5  is  a  conventional  differential  expression  representing  one  degree 
of  freedom.  On  the  right,  we  have  the  force  in  the  system  immediately  above^it  (weighted 
by  the  mass  ratio  ),  and  an  acceleration  forcing  term  equal  to  the  input  ,  less  the 
accelerations  of  the  springs  between  the  point  of  application  of  and  the  system  under 
consideration. 

For  convenience  we  can  write 


(7) 


so  that  equation  5  becomes 


I 


(8) 


We  can  rearrange  these  equations  by  substituting  into  the  right  side  from  the  previous 
equation,  so  that  we  ejqjress  the  motion  of  a  particular  system  in  terms  only  of  the  two 
systems  on  each  side  of  it. 

Thus, 


(9) 

Sf»  ^.1 )  ^  *  ^-1 


The  relationships  given  in  equation  9  are  of  particular  importance  in  analog  hookups, 
where  we  wish  to  minimize  circuit  complexity.  The  basic  flqw  chart  is  shown  in  figure 

3. 


b.  The  effect  of  parallel  elements. 


If  some  elements  occur  in  parallel,  the  generalized  equations  become  coupled  in  the  fol¬ 
lowing  manner,  using  the  notation  of  figure  4. 


t*ti 

h 


5. 

% 


Figure  4.  A  branched  system. 


Below  the  branch 


0  F*-  =  ^  Ff,  F^^ 


Above  the  branch,  in  the  system 

%  *  ''h  ’  ^ 

•*  j  ^ 

The  equations  are  the  same  in  the  ^  branch,  ^  being  substituted  for  ja  . 

c.  General  equations  for  a  series  degree  of  freedom  without  mass. 

A  "restraint  system"  can  often  be  represented  as  a  dynamic  system  in  series  with  a 
dynamic  model.  The  mass  of  the  restraint  system  is  considered  negligible. 

If  the  restraint  system  mass  is  finite,  it  can  be  regarded  simply  as  system  (1)  in  figures 
2  or  3.  To  study  this  case,  consider  the  first  expression  of  equation  5: 


i,  ^ ^ ■i’ 

by  making  the  substitution 

Then 


Now  by  definition, 

Force  in  damp';r  0)  =  f.a)  -  ■p. 
and  force  in  spring  (t)  ~  ^ 


i 


(It  should  be  noted  that  and  are  actually  functions  of  and  1**,  respectively, 

so  that  and  Sy  are  independent  of  . ) 

Substituting  into  (12),  we  have 

Now  putting  =  0  gives  ^  -f  »  “Px  -f 

'  f,  = 

We  can  approximate  this  condition  on  an  analog  by  making  large  compared  to  unity; 
of  the  order  cf  100  sa3^  Then  the  frequency 

JL  -  A'/jl) 

^  fC&/ 

A  linear  restraint  system  can  "bottom"  however,  and  subsequent  to  this  the  acceleration 
input  is  fed  directly  to  the  dynamic  model,  together  with  an  impulsive  velocity 

change  which  is  equal  to  the  velocity  at  which  the  restraint  spring  bottomed.  This  latter 
requirement  involves  separation  and  the  differentiation  of  velocity  signals,  and  it  is  found 
that  for  analog  work,  more  conventional  methods  of  programming  are  preferred. 

For  nonlinear  restraint  systems,  however,  this  system  can  be  used  very  effectively  since 
there  is  no  discontinuity  in  the  equations  of  motion. 

When  the  mass  of  the  restraint  system  can  be  neglected,  the  dynamic  model  is  as  shown 
in  figure  5.  Equation  8  still  applies  to  this  system,  except  that  the  terms  for  n  1 
are  replaced  by  the  simple  force  equality 

f,  =Fi  (13) 

Thus  the  equations  become 

F,  ^  Fx 

C ye 


(14) 


S,  is  an  unknown  which  must  be  determined  from  the  first  equation.  Substituting  equation 
7  for  and  ^  , 


»  id,) 


By  differentiating  this 
for  $f  in  terms  of 


equation  with  respect  to  time  it  is  often  possible  to  obtain  a  solution 


Figure  5.  Zero  mass  restraint  system  in  series 
with  a  dynamic  model. 


From  the  foregoing  it  is  clear  that  we  may  express  the  influence  of  zero  mass  restraint 
system  in  the  following  theorem. 

"The  only  influence  of  a  zero  mass  restraint  system  on 
the  equations  of  a  dynamic  model  is  to  replace  the  forcing 
term  Ij^  with  the  term  )>  wherever  it  appears.  " 


i 


2. 


:  INEAR  SYSTEMS 


For  a  linear  system,  let 


- 


damper  force  ^ 


(t,) 


C,  -  i>  -  ^sli 

Mf 


kr 


spring  force 


k 


CO'  *•  Jlf  ss. 

^4- 


(S^) 


Then  equation  5  becomes 


^/v/  ^y*+/ } 


Equation  9  becomes 


%  0^  L,  -f-  ?U,] 


(16) 


(17) 


m  • 

When  the  excitation  (  ^  )  is  sinusoidal,  or  can  be  resolved  into  a  Fourier  series, 
equation  17  can  be  solved  using  conventional  impedance  techniques.  Alternatively, 
we  can  use  Laplace  and  Fourier  transformations,  v/hich  convert  the  differential 
equations  into  algebraic  equations.  For  example  the  Laplace  transformation  of 
equation  16  is 

^  r  c  — ' 

-i-  initial  conditions.  (18) 


10 


This  set  of  simultaneous  equations  can  be  solved  for^  ). 
for/  ( $ )  can  be  inverted  to  S  =  -f  ( ^ )  by  means  of  tables 
the  integral 


An  analytical  expression 
of  Laplace  transforms ,  or 


where 


denotes  integration  about  a  closed  contour. 


11 


3. 


CLOSED  FORM  SOLUTIONS  FOR  SPECIAL  LINEAR  CASES 


a.  V.Tiolly  viscous  restraint. 


In  this  section  we  consider  the  idealized  case  of  a  restraint  system  whose  mass  and  spring 
characteristics  are  negligible  compared  with  its  damping  characteristics,  as  indicated  in 
figure  6.  I — - , 


Figure  6.  Wholly  viscous  restraint  system. 


Equation  16  then  simplifies  to 

# 

S»  * 


(i<  U) 


(19) 


Differentiating  the  restraint  equation 


c,  *  2«;  * 

•  • 

Substituting  for  in  the  first  c  f  equation  19,  the  equation  of  motion  becomes 


5, 


(/  y-<i/<r,) 


<0. 


= 


(20) 


(21) 


We  have  thus  reduced  a  two  degree  of  freedom  problem  to  a  single  degree  of  freedom 
equation  in  S-  .  We  know  that 


+ 


C%/C, 


(22) 


the  motion  will  be  subcritically  damped. 


2 


Figure  7.  Restraint  sybtijm  clamping  required  to  give 

"Dead-Beat"  motion  of  linear  djmamic  model. 
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Insofar  as  dynamic  overshoot  of  the  human  body  is  a  problem,  at  least  for  short  period, 
short  rise  time  pulses,  a  ’’deadbeat"  restraint  system  is  of  considerable  interest.  From 
equation  22  this  could  be  achieved  if 


where  *  CA**  J  ®  A;* 

gx 


(23) 


Eqpiation  23  is  plotted  in  figure  7,  and  a  table  of  values  is  given  below. 


c 

.0 

.25 

.  6 

.  156 

.05 

.238 

.7 

.147 

.1 

.228 

.8 

.139 

.2 

.208 

.9 

.132 

.3 

.192 

.  95 

.128 

.4 

.179 

1.0 

.125 

.5 

.167 

Closed  form  s®lutions  to  the  complete  problem  cannot  be  obtained  because  it  is  impossible 
to  obtain  a  closed  form  solution  to  the  equation  for  time  to  bottom-out  the  restraint  system. 
Th^s  point  is  best  illustrated  by  considering  the  even  simpler  system  of  figure  8,  where  the 
damping  of  the  human  body  dynamic  model  is  assumed  to  be  zero. 


This  simplification  peimits  the  auxiliary  (restraint)  equation  to  be  written  as 

>(p,  *  ttC, 


(24) 


Figure  8. 


Dynamic  system  for  a  zero  damping  model. 


Defining  /***%.  ^ 


tAy 


or 


t,-  #/ 


c^, 

eti 

t 

ix 


The  basic  equations  of  motion  of  the  system  are  then 

S,  V  37  S.  (*<U) 


(25) 


(26) 


(27) 


where  “tg  is  the  time  at  which  the  restraint  system  bottoms  out  (  ^^  * 
In  terms  of  a  nondimensional  time  't  to’t' 


,  ys 

dJt 


The  equations  of  motion  become 


At'- 


^r‘ 

+  =  f 

A 

Crcr,) 

+  $X  ^ 

and  from  equation  26 

i,  = 

3-7  /  Jt 

(28) 


(29) 


From  figure  7  we  see  that  if 


C  <  1/4,  the  damping  is  subcritical. 
g  ^  1/4,  the  damping  is  supercritical. 


(1)  Response  to  an  impulsive  velocity  change  . 


When  the  forcing  function  is  zero,  and  the  initial  velocity  (  5 


has  the  solutioi. 

t'ir 

^  s;«,'vr 

7 

^M4ft 

s  e 

)  =  ,  equation  28 

4 


(30) 

(31) 


..  -/ 


where  ^  ^  in  the  second  quadrant. 

Also,  from  equations  29  and  30 


(32) 


We  obtain  the  maximum  value  of  S,  by  differentiating  equation  32  with  respect  to  . 

-rA^ 

s')  (33) 

where  Pa  +  ‘•j 


1  e, 


or 


ATT 


'f 


5,  is  a  maximum  when  f' 

V  -  itp- 


] 


(34) 


It  13  therefore  possible  to  plot  curves  of  ^  J a  function  of  C.  ,  and  also 

curves  of  rom  these  we  can  determine  the  variation  of 

with  allowable  bottoming  depth 


1  e 


co/i>o~ 


eA> 


'ha* 


-  it/'- 


-2 

e 


7 


(35) 


(36) 


(2)  Response  to  a  zero  rise  time  acceleration 


•  • 

yi  . 


For  this  case  the  right  side  of  equation  28  becomes 

y.  = 


The  solutions  are 


-v/^r 


X*-  -  /  -  — 

t  '  ^ 


-f-f) 


(37) 


where 


f 


•  -  r 
«  5<t* 


(38) 


The  restraint  deflection  is 


i' '  tcli'-  f 


dX 


and  after  some  manipulation  this  becomes 


(39) 
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In  this  case  there  is  no  true  maximum  deflection  since  S;  continually  increases,  and 
it  is  therefore  logical  to  define  the  critical  bottoming  depth  as  the  value  of  when  S 
reaches  its  maxinium. 


That  is,  when 

^  f  * 


where  ^  *  *"  ^  T 


t  • 


4c? 


htt 
^^*7 . 


c 

«  •  ^ 


«  o 


(40) 


(41) 


Substituting  equation  40  for  “Xm.  in  equation  39 


9< 


M/'ir  ^ 

*7  "  2?/ 


e 


(42) 


Note  than  when  C  »  1.0  7->l,0 

c 

^  2.0 

b.  Linear  Spring  Restraint. 

The  corollary  to  the  idealization  of  Section  a.  is  the  system  illustrated  in  figure  9. 
For  this  case  equation  16  becomes 


^  4-  -  jc  - 

I,  +  Si.  * 


(43) 
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(a) 

O) 


Figure  9.  Linear  spring  restraint. 


Differentiating  the  restraint  equation 


s, 


+■ 


(44) 


The  presence  of  the  third  derivative  makes  any  solution  difficult,  so  we  will  once  again 
use  the  zero-damping  model,  by  writing  ^  =  0. 


Defining 


and 


(1)  Solution  for  shoi*t  period  ( ^  t  >  at,  )  acceleration  with  zero  rise  time. 

•  • 

We  assume  that  the  acceleration  inrnit  rises  instantly  from  zero  to  *  Y.  Then  from 
equation  I.  3  in  Appendix  I 


^  / 

A/tTf 


/  - 


(47) 


(48) 
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let  -  f 


Now 


t 


«Oz 


P  w 


/J7T^  /i/l-t  ^  A^l-i-^ 


J> 


Si 


g 

•  • 


^f-tf 


/UlTf 


(49) 


and  of  course,  from  equation  46 


Wf 


so  that 


'"jih*  =  J2.t  -  f  ‘ 

y: 


(50) 


Therefore  the  total  (initial  condition)  spring  velocity  after  bottoming  will  be 


Y*  ^ 


(51) 


Substituting  these  initial  conditions  into  equation  1. 6,  in  Appendix  I,  the  motion  after 
bottoming  is  described  by  the  equation 


«  t  —Jt*  t  i>) 


y< 


where  S 


M-U  ^  « 


-U-  S} 


//t-h  fltt-il 


(52) 


(53) 
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Obviously  the  maximum  value  of 


» r 

•  # 

Y, 


”  •  • 

Y. 


IS 


=  /  ■*■  jJi+fiU-  fJ 


An  interesting  limit  is  provided  by  the  case  0. 

c 

Here  <0|t 


/  +/i  * 


This  checks  with  other  zero-zero  cushion  solutions. 


(54) 


(rib) 


For  convenience  let 


/  -  -  m':) fi 


(56) 


Then  f 


-  (  -f-  4/ 1  -f  'f£i-  'f/fl 


The  critical  value  is 


(57) 


*  2 


(58) 


For  y  >  ■V'«rr  equation  57  is  not  applicable,  and  the  solution  is  merely 

^  r 

2  d« 


2.0 


P  p 

K 


The  maximum  value  of  equation  57  is  found  by  differentiating  with  respect  to  and 

equating  to  zero,  holding  £l>^  constant. 


1  e 


/  /yt  ^ 

(  /-  'TfT/ 


X  =  -4i- _ 

y,  /  V t+i'U-Wf^ 


•=  O 


'  f 


(51)) 


Substituting  into  equation  L  2,  Appendix  I 


-  / 


( 

; 


Also,  when  the  restraint  stiffness  is  zero 


f 


and 


^  ^  ^  a/ f 

I  e 

This  is  the  "slack  harness"  (or  zero  cushion  stiffness)  solution. 
(2)  Impulsive  velocity  change  ^  ^  . 


(60) 


(61) 


For  an  impulsive  velocity  cliange,  we  can  regard  the  problem  as  a  single  degree  of  free¬ 
dom.  Prior  to  bottoming  the  effective  spring  stiffness  is 


X  * 


#  * 


From  Appei 


I, 


Jl 


(62) 

(63) 


(64) 


For  equality  of  force  in  the  springs, 
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c  » 

and 

^  Sj- 
/■h  f 

(66) 

Thus ,  from  equations  64  and  66 

= 

4^.  slu.  SLt. 

SL  * 

(67) 

JL  ^  ^ 

(68) 

• 

(69) 

i  ^  f 

(70) 

After  bottoming,  the  spring  stiffness  is  kt.  and  the  velocity  of  the  mass  with  respect  to 
the  ground  is  Sj.  .  Thus,  equation  64  becomes 


where  £<  =  the  energy  storage  parameter,  being  the  total  energy 

the  resiliency  can  absorb  at  bottoming  (ft.  lb. ) 


^  =  mass  of  occupant  (slugs) 


The  most  favorable  condition  (  a  minimum)  is  for  the  restraint  system  to  just  bottom 
out  under  the  influence  of  an  impulsive  acceleration.  This  occurs  when 

Aat  a 

or 


c.  "Crushable  foam"  (Constant  force)  restraint. 

This  linear  case  may  likewise  be  regarded  as  the  simplest  example  of  a  nonlinear  restraint 
system,  and  an  explicit  solution  obtained. 

We  divide  the  system's  behavior  into  three  distinct  regimes: 

(A)  when  the  spring  force  kt  foamds  "rigid.  " 

(B)  whenyf'j^fjs  F  the  foam  deflects  but  the  spring  does  not. 

(C)  when  the  foam  bottoms  out  it  becomes  rigid  again,  and  in 
addition  to  the  acceleration  input  ,  the  dynamic ^model 
has  to  accommodate  an  impulsive  velocity  change  • 

9  » 

We  now  evaluate  these  regimes  for  a  short  period  acceleratio  i  of  magnitude  and  zero 
rise  time. 


Figure  10. 


Dynamic  model  for  a  constant  force 
restraint  system. 


(1)  Short  period  acceleration  solution. 


Condition  (A) 
For  this  regime 


c  +  -  r. 


(76) 


If  the  initial  spring  deflection  is  zero,  the  solution  to  this  equation  is 

=  t  - 

Obviously,  the  end  of  this  regime  is  defined  by 


y. 

9 

^  * 

/—  ^Sht^  tO* 

or  "i 

SS  /  — 

The  spring  velocity 

at  this  time  is  given  by 

(77) 


(78) 


(79) 


Condition  (B) 

The  spring  has  now  reached  its  critical  deflection  and  remains  at  this  value. 

S,  = 

Si  =  At  -  j/) 


(80) 


But 


‘ji  ^  F/m 
. 

•  • 
S, 


and  =  0 

y,  +  L  ' 

••  ••  «  M  r 
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(81) 


c  -  ^  y  "1  -4  ^ 


Bottoming  out  occurs  when  S(  -S,t 


1  e, 


- 


thsf  Ft  -  g 
%-Ff-^  y.-F/^ 


5%  \-'^J'^  T"Cy^- 

X-  F/^  I  ^  ^tCFL  ^ 


Substituting  equation  79  for  and  defining 


* 


gMA 


X 

f_,  ^ 


(82) 


(83) 


(84) 


Note  that  from  equation  78 


C(rt 


4..  = 


5«Ca. 


(85) 


/-/ 


(86) 


From  equation  81  the  foam  will  bottom  out  with  a  velocity 


«  * 


fCg-P 


(87) 


Condition  (C) 


In  condition  (C)  the  equation  of  motion  returns  to  the  form  of  equation  76.  The  initial 
conditions  are  no  longer  zero,  and  the  solution  becomes 


Now  Sz  =  0  prior  to  bottoming,  and  the  velocity  change  Sa  has  been  acquired  im¬ 
pulsively  by  the  spring  at  the  instant  of  bottoming,  thus  substituting  equation  86  into  88, 
and  remembering  that 

^  -§)  *}0~ 

(90) 


^  id-i) 

'-i  0-iT 


(91) 


where  is  given  by  equation  86. 

If  we  define  ^  ■*“  ^ ^ f  “  ^') 


(92) 


(93) 
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Thus,  the  only  variables  are 


and 


I 


Std 


Some  limit  cases  are  of  interest.  For  example,  as  ^  *  "y 


OtMA* 

9  9 


1.0, 


2>0 


Similarly  for  ^  0 


LS 


I  -f  //* 


•  • 

which  is  the  same  as  eciuation  55. 

(2)  Impulsive  velocity  change^ 

i 

In  this  case  the  energy  to  be  absorbed  by  the  total  system  is  .  Before  the  foam 

starts  to  crush,  the  model  will  absorb  ^ 


^^0 


(94) 


If  ^  is  the  crushing  force,  and  the  bottoming  depth,  the  energy  absorbed  in  the 

foam  will  be  F  .  Thus,  the  remaining  energy  to  be  absorbed  by  the  spring  will  be 


z 

0 


(95) 


Dividing  throughout  by  M  , 


(96) 

(97) 


^  9m 

Naturally,  if  the  initial  energy  2  is  insufficient  to  bottom  the  foam,  then 


Also,  if 


< 


(98) 


The  result  will  be 


4J2 


(99) 


A  convenient  parameter  is  therefore 


(100) 


The  optimum  crushing  force  for  a  given  velocity  change  is  the  one  which  results  in  the 
lowest  value  of  for  a  given  bottoming  depth  • 

This  condition  obviously  occurs  when  the  foam  bottoms  at  the  same  instant  that  all  the 
initial  kinetic  energy  is  absorbed.  In  other  words,  from  equation  95 


^  M  — - 

f  ~  i 

(101) 

But  ^  ^0 

- 

^  /ft 

(102) 

/.  ^z*’  ^0) 

(103) 

and  from  equations  101  and  102 

(104) 


and  Dynamic  Response  Index  (DRI) 


UKl  with  optimum  restraint  _  —  ^  ! f  ^  I 

DRI  with  no  restraint  A'  ^  '  (105) 


d.  The  influence  of  slack  in  the  restraint  system. 


When  there  is  slack  in  a  restraint  system,  or  the  occupant  is  separated  from  it,  as  in 
ejection  from  an  aircraft  subject  to  negative  g,  an  initial  velocity 


(106) 


is  built  up  prior  to  contact  with  the  restraint  system  (which  occurs  at  time 't'  =  jf"* ). 
Thus,  the  equations  already  developed  are  still  valid,  but  the  initial  velocity  condition 
also  appears  in  the  solution,  increasing  (most  usually)  the  physiological  effect  of  the 
applied  acceleration^ 


The  "slack  distance" 


is  related  to  the  input  acceleration  by  the  equation 


(107) 


In  analog  work  Sf  can  be  deterr  aned  as  a  function  of  for  any  input  acceleration- 
time  history.  The  value  of  appropriate  to  the  problem  likewise  can  be  determined 

and  fed  into  the  computer  as  an  initial  condition. 


It  is  clear  that  if  the  acceleration  input  is  an  Impulsive  velocity  change,  restraint  system 
slack  will  have  no  influence. 


(1)  The  effect  of  slack  when  the  restraint  system  is  rigid. 


For  a  rectangular  (zero  rise  time)  acceleration  pulse  (  )  in  the  short  period  regime, 

the  velocity  built  up  before  contacting  the  rigid  restraint  system  will  be 


(108) 


and  the  travel 


i 


(109) 
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Thus 


(110) 


t, » yz 


and 


(111) 


Substituting  this  initial  condition  in  equation  L  7 ,  Appendix  I 


S 

T 


Tf  = 


I  £-  Jt-  +  -f  z) 

7 


where 


(112) 

(113) 


(114) 


We  now  have  to  find  the  solution  for  the  maximum  value  of  equation  113  which  can  be 
done  using  the  techniques  of  section  L  3  in  Appendix  1.  A  much  simpler  solution  can 
be  obtained  for  zero  damping,  of  course,  since  equations  113  and  114  then  become 


(115) 


(116) 


is  a  limit  solution  and  has  been  obtained  elsewhere  in  this  report.  Equation  116  is 
plotted  in  figure  11.  It  is  interesting  to  note  that  a  slack  of  only  half  an  inch  (  =  .  04167) 

in  the  spinal  mode  ( =  251. 0  rads/sec)  would  increase  the  DRI  of  a  20  g  acceleration 
pulse  (with  zero  rise  time)  by  as  much  as  100%. 


(2)  The  effect  of  slack  in  a  linear  spring  restraint  system. 

The  equation  of  motion  for  this  case  is  given  in  Section  3,b.  When  damping  in  the  dynamic 
model  is  neglected  equation  45  gives  the  solution 


(117) 


Figure  11.  Increase  in  D.  R.  I,  due  to  slack  when  the 
restraint  system  is  rigid,  for  a  zero  rise 
time  pulse.  (Zero  damping  dynamic  model). 
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After  the  spring  has  bottomed  we  have,  from  equation  f)! 

•  •  X 

X 

Fi'om  equation  117 

«z£„  ,  je*.  si*,  ^  /B  ^ 

yTff  A^f 


(119) 


(120) 


From  equation  118 

^  f -il 

Obviously  but  the  different  signs  in  equations  120  and  121  prevent  us  from 

substituting  for (&  i- ^ ) in  equation  120.  Thus,  a  single  explicit  relationship 
cannot  be  obtained. 

By  substituting  eouation  120  into  119  we  have 


a  a 


(122) 


and  is  given  by  equation  121. 

e.  The  influ^^ace  of  preloading  in  a  restraint  system. 


The  effects  of  preloading  can  be  examined  effectively  with  an  analog  computer,  where 
preloading  can  be  simulated  as  an  initial  condition  acceleration  input 


For  the  simple  case  of  a  rigid  restraint  system,  however,  subjected  to  a  zero  rise  time 
acceleratior  we  can  obtain  a  solution  from  the  work  in  Appendix  I.  From  equation  1. 7 , 


f 


=  0,  so  that  Thus,  equation  1. 14  becomes 


%  • 


'<»<»»■  =s  t  -1-  f 


iTc.  t  C  V 


(123) 


In  figure  12  the  dynamic  overshoot  for  zero  preload  is  plotted,  and  in  figure  13  the  in¬ 
fluence  of  preload  is  included.  It  was  found  in  the  last  section  that  the  parameter  in¬ 
fluencing  the  physiological  effect  of  a  zero-rise  time  acceleration  was 

where  $$  is  the  slack  distance.  For  preloading  the  parameter  is 


(T 

*'2* 


y 


preload  force  in  lb,  resolved  along  the  appropriate  axis 
weight  of  occupant  in  lb  x  acceleration  in  g's 


As  an  example^  consider  the  spinal  mode  which  is  of  importance  in  aircraft  ejection  seat 
work  and  for  which  C*  =  0. 224.  The  maximum  tolerable  preload,  even  with  a  powered 
inertia  ree?,  is  probably  about  1  g,  resolved  along  the  spinal  axis,  so  that 


5'Oi:' 


«  • 


=  15  g  say, 


K 


.0667 


Thus,  clw  attenuation  factor,  from  figure  13,  is  about  0.98,  or  a  2%  reduction  in  DRI. 
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WITHOUT  PRE-LOADING 


SESTHAINT  PRE-LOAD  PARAMETER  » 

^ _  PRE-LOAD  FORCE  (LB) 

IT  WT.  OF  OCCUPANT  (LB)  X  ACCELERATION  (^*s) 

Figure  13.  Influence  of  preloading  restraint  system  on  of 
(Damped  linear  system  ^bjected  to  a  zero  rise 
time  acceleration  input  ^  ). 
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It  is  obvious  from  the  foregoing  that  preloading  the  restraint  system  has  little  value, 
from  a  physiological  point  of  view.  Since  it  is  of  extreme  importance  to  avoid  slack, 
however,  as  shown  in  the  preceeding  section,  some  preloading  should  always  be  in¬ 
troduced  in  the  process  of  cinching  down  the  restraint  system  to  make  sure  that  no 
slack  exists. 
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4. 


SOME  SIMPLE  NONLINEAR  RESTRAINT  CASES  INVOLVING  CONTINUOUS 
FUNCTIONS 


The  simplest  type  of  nonlinear  restraint  problem  consists  of  a  zero-damping  model  of 
the  human  body  in  series  with  an  undamped,  nonlinear  spring.  Thus ,  the  problem  may 
be  regarded  as  a  nonlinear  extension  of  the  work  in  Section  n,  3,  b. ,  page  18. 


From  equation  14  the  appropriate  equations  of  motion  will  be 


or,  more  specifically 


(124) 


(125) 


where  ^  (  /,  )  is  some  continuous^function.  By  differentiating  the  {  Sf)  equation 
we  can  obtain  an  ejq)ression  for  which  in  some  (very  few)  cases  enables  a  general 
solution  to  be  obtained. 


Retaining  the  generality  of  equation  125 


.  c  = 

^$1  ^ 


Thus, 


(126) 


(127) 


(128) 


Substituting  for 
,  by  a  prime? 


,  from  equation  126  and  denoting  differentiation  with  respect  to 

_  f(a) 


(129) 
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Figure  14.  Simple  model  of  a  nonlinear  restraint  system. 


Substituting  into  equation  1.25 


(130) 


Even  though  we  have  succeeded  in  reducing  the  problem  to  a  single  degree  of  freedom, 
this  equation  is  difficult  to  solve,  whatever  the  nonlinear  function  (  S,  )  niay  be. 
The  only  cases  which  can  be  solved  this  way  are  of  little  value  in  practical  work. 


As  an  alternative  approach,  we  can  consider  the  two  springs  in  figure  14  as  being  one 
spring  with  a  rather  more  complex  force-deflection  relationship.  This  approach  is 
more  promising  because  the  solution  for  a  single  nonlinear  spring-mass  problem  can 
be  fairly  simple. 


Retaining  our  generality,  and  defining  the  total  deflection  of  the  two  springs  as 


(131) 


If  this  can  be  solved  for  in  terms  of  Sir  ,  then  the  resulting  equation  describes 

the  equivalent  single  nonlinear  spring. 


For  the  general  nonlinear  restraint  spring  examined  in  the  next  section,  equation  131 
takes  the  following  form 


JL 

^  'A  (  014/ 


(132) 


(133) 


and 


2. 


^7 


3i 


(134) 


In  order  to  evolve  an  explicit  expression,  we  must  solve  this  polynominal  in  (  FA), 
which  is  easy  to  do  for  ivt  =■  1  or  2.  For  =  3  or  4  the  solutions  become  increasingly 
complicated,  and  higher  powers  become  impractical. 


A  closed  form  solution  is  possible  for  the  simplest  polynominal  nonlinearity,  which  is 


£■  = 

r 


t 


(135) 


Since 


^2 

£.  «  St  *“^0  ^  { Sr 


(136) 


and  the  solution  for  (  /%<)  is  found  to  be 


F 

tfM 


(137) 


The  sane  type  of  solution  could  also  be  obtained  for  the  cxibic  nonlinearity 


(138) 


but  with  added  complexity. 

Another  type  of  nonlinearity  considered  in  the  next  section  is  the  tangent  relationship 


(139) 


(140) 
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Obviously  no  explicit  expression  for  can  be  obtained  in  this  case. 

From  the  foregoing  it  is  evident  that,  even  for  the  simplest  possible  zero  damping 
cases,  it  is  generally  not  possible  to  achieve  closed  form  solutions  for  the  effects  of 
a  nonlinear  restraint  system,  except  for  the  quadratic  case.  Consequently,  the  only 
practical  method  of  treatment  is  to  employ  an  analog  or  digital  computer. 

a.  Generalized  power  law  nonlinearity. 

When  the  spring  force-deflection  curve  is  described  by  the  equation 

the  equation  of  motion  of  a  simple  system  becomes 

(1)  Solution  for  an  impulsive  velocity  change 

The  potential  energy  stored  in  the  spring  at  a  deflection  SmAx  is 

^  r -  s.  1 


(141) 


(142) 


(143) 


(144) 


,  .  t 

Equating  144  to  the  initial  kinetic  energy  g  m  #Ur  before  impact, 

0\+f  L  J 


=  f  r  ^  n 


(145) 


For  zero  initial  conditions,  the  DRI  is 


7“  =  7-  (  t) 


(146) 
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(2)  Solution  for  constant  short-period  acceleration 

# 

Let  $  -  in  equation  142 

r  as 


C  ^  h 

X  liCU  ff  ^  - 


(147) 


For  zero  initial  conditions 


^  .  ii. 


and  the  DRI  is 


(148) 


b.  Generalized  cubic  non-linearity. 

In  this  section  we  are  concerned  with  an  undamped  model  whose  spring  force-deflection 
curve  is  described  by  the  equation 


(149) 


so  that  the  equation  of  motion  becomes 

$  ^{JS 


(1)  Solution  for  an  impulsive  velocity  change  . 

The  potential  energy  stored  in  the  spring  at  a  deflection  is 


£ 

£_ 


9 


(150) 


(151) 


(152) 
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(153) 


Equating  152  to  the  initial  kinetic  energy  before  impact 

s  -  SI*)  ^ *"  ) 

For  zero  initial  conditions 


The  DRI  is  obtained  by  substituting  for  i*'  equation  149. 

(2)  Solution  for  constant  short-period  acceleration  Vc  • 

m 

Let  S  ~  in  equation  150 
Then 


4P. 


1  (L..  -s.)-  tVc.  -s:)  -i  (C.  -s:) = 


(156) 


For  zero  initial  conditions 


•  • 


At  A* 


-ifS. 


MAt 


O 


(157) 


O  is  given  by  the  solution  to  this  cubic. 

ptA* 

c.  Generalized  tangent  law  nonlinearity. 


■n _ _  J-I.J..-  ^  ^  #%  tir»4V*  ri  o  fl-k-rvi  v«rr  r»^' r*  1 1  n  OO  T*  (r*V»f!l  T'Q  r*"#  fiT'l  g1  1  r'  IG  G 

Wl?  ill  c;  VTXVii  1*  o|./A  T»»»VWUW  >«v  .  v  ^  v  * 


defined  in  figure  15. 
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Figure  15. 


Tangent  spring  force  deflection  curve. 


1  e 


F-  tfc:;  211^ 

-tr  2 


(158) 


Writing 


f  -  - 


(159) 


the  equation  of  motion  is 

g  ^  jb  fs  =  y. 

mg 

k  ^ 

or  writing  — 


(160) 


s  ^  ff  =  y 

£ 

(1)  Solution  for  an  impulsive  velocity  change 


The  potential  energy  stored  in  the  spring  at  a  deflection  is 

*  c  •  ^ 


7*  Sc  ! 

Equating  162  to  the  initial  kinetic  energy  ^  ^  J/tr*  before  impact 

-I--  -  -f-'it/sfrV 


(161) 


(162) 


(163) 


for  zero  initial  conditions 


(7  '  ' 


^  f 


^<4*' 


(164) 


and  the  DRI  is 


1[ 


MA* 


Now 


K  = 


r  -  ( 


F4 


MA*  ^ 


AA 


(2)  Solution  for  constant  short  period  acceleration  Vc 

« 

Let  5  =  />  in  equation  161 


Then 


p  0 

•'■  J '//I '  fluffs)! 

•  •'iS.  ?  ^ 


(166) 


(166) 


(167) 


(16b) 


For  zero  initial  conditions 


-6*>- 


P<A!t 


(169) 


There  is  no  explicit  solution  for  the  DRI  in  this  case. 
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SECTION  III 


THE  DERIVATION  OF  DYNAMIC  MODELS  OF  THE  HUMAN  BODY 
1.  GENERAL  OBSERVATIONS 


The  basic  aim  in  Body  Dy  lamics  is  to  devlve  a  system  of  differential  equations  which 
are  in  harmony  with  the  observed  behavior  of  the  human  body  when  it  is  subjected  to 
acceleration.  Since  there  is  always,  or  -iearly  always,  a  physical  analog  to  a  differ¬ 
ential  equation,  it  is  frequently  easier  to  discuss  the  analog  rather  than  the  equations 
it  illustrates  because  our  thinking  is  typically  better  adapted  to  considering  physical 
objects  than  to  reasoning  in  terms  of  mathematical  abstractions. 


A  typical  "dynamic  model"  is  illustrated  in  figure  16.  This  analog  has  two  degrees 
of  freedom  because  the  masses  can  move  in  an  arbitrary  manner  with  respect  to  each 
other. 


INPUT  ACCELERATION 


Figure  16.  Dynamic  analog  or  model  for  a  two 
degree  of  freedom  problem. 


The  model  used  in  figure  16  is  a  "lumped  parameter"  model.  That  is,  the  mass,  spring 
and  damping  characteristics  are  all  considered  as  separate  and  discrete  elements.  Such 
treatment  is  not  always  possible.  Dynamic  models  of  the  head,  developed  to  predict  head 
or  brain  injury,  are  based  on  the  assumption  that  mass,  damping  and  stiffness  are  distrib¬ 
uted  uniformly  throughout  the  skull.  As  might  be  expected,  the  equal  ons  for  this  case  are 
rather  more  complicated  than  for  an  equivalent  lumped  parameter  model,  but  tL^^re  is  al¬ 
ways  a  lumped  parameter  equivalent  of  such  a  system. 

Dynamic  models  of  the  human  body  are  based  mainly  upon  experimental  measurements  of 
its  response  to  various  types  of  acceleration  inputs  as  well  as  the  subjective  reactions  of 
the  individuals  to  the  tests.  Experiments  with  cadavers  are  also  of  value,  although  to  a 
lesser  extent  since  substantial  changes  occur  in  some  of  the  body's  dynamic  characteristics 
after  death. 


Perhaps  the  most  cogent  reaso»:  for  the  ust  .jI  dynamic  models  is  thrr  I'ley  >;royide  an 
essentially  unifying  theory  to  physical  u)'»scr'n\tions  of  human  resiJOiu  -:  u  aceeler;.t;  -.!: 
thereby,  permitting  impaci  data,  nivcraft  ejection  seat  firings,  sled  tests,  siiipboard  '/n  ■ 
juries,  and  a  host  of  other  observations  to  he  reduced  to  a  common  /dr  ator,  kmwn 
as  the  "Dynamic  Response  Index"  (DRI).  A  theory  ftliich  enables  all  experimental  c  h. 
vations  to  be  reduced  in  tenns  of  a  single  parameter  can  also  be  used  in.  reverse,  of 
course,  to  determine  the  probable  results  of  an  experiment  which  has  yet  to  be  carries! 
out;  or,  putting  it  another  way,  to  determine  whether  the  experiment  is  worth  carrying 
out  at  all. 

In  the  present  context,  we  are  concerned  with  determining  the  numerical  constants  of 
lumped  parameter  dynamic  models  which  will  best  represent  the  human  body.  We  use 
as  the  basis  of  our  effort  the  available  experimental  data. 


Experimental  data  can  be  divided  into  three  classes; 

(a)  Experiments  in  which  the  acceleration  input  pulse  is  of 
a  type  which  enables  the  data  to  be  used  directly  in  de¬ 
termining  the  dynamic  model  parameters ,  The  best 
example  in  this  class  is  an  impulsive  or  "Impact"  velocity 
change. 

(b)  Cases  in  which  the  data  must  be  reduced  by  means  of  a 
dynamic  model  before  it  can  be  used.  An  irregular 
ejection  seat  acceleration-time  history  is  an  example 
of  this  second  class. 

(c)  Cases  for  which  some  or  all  of  the  acceleration  records 
ar''  \v.  .thlpsR.  but  for  whic*"  an  pppropriate  idealisation 
of  the  acceleration  input  can  be  deduced  indirectly. 

In  general,  it  is  easiest  to  use  type  (a)  data  in  the  generation  of  a  dynamic  model,  and 
then  to  use  this  model  to  improve  probability  of  injury  estimates  with  type  (b)  and  (c) 
data. 

There  are  four  different  types  of  information  which  are  of  value  in  determining  the  values 
of  a  dynamic  model  for  any  particular  direction. 

From  acceleration  experiments  we  can  determine: 

(a)  the  critical  impact  velocity  change  , 

•  • 

(b)  the  critical  short  period  acceleration  .  and 

(c)  the  "free  oscillation"  frequency  (if  the  subject 
is  unrestrained  in  I’ebound). 


From  impedance  measuremento  v/e  can  determine; 


(d)  the  "small  amplitude"  resonant  frequency 
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2. 


METHODS  OF  DETERMINING  LUMPED  PARAMETER  MODELS  FROM  THE 
AVAILABLE  DATA. 


a.  Analysis  of  acceleration  tolerance  data  alone. 


Experiments  in  which  an  approximately  rectangular  acceleration  pulse  (figure  17)  is 
imposed  on  a  live  human  subject  should  result  in  iujur;  (or  more  precisely,  jirobability 
of  injury)  curves  of  the  type  shown  in  figure  18. 
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Figure  17.  Definition  of  rectangular 
acceleration  pulse. 


Figure  18.  Typical  injury  curves  for 
rectangular  acceleration 
pulse. 


The  "corner  duration"  )  marks  the  boundary  between  two  dynamic  regimes.  For 
are  in  the  impact  or  impulsive  wl  'city  ch'^nge  regime.  The  pulse  shape  has 
no  influence  upon  the  behavior  of  the  system,  but  only  the  pulse  area,  which  is  equa.1  to 
the  velocity  change  imposed.  In  practical  terms,  the  pulse  duration  in  this  regime  is  so 
small  that  it  is  over  before  the  body  has  started  to  respond.  Thus,  the  body  sees  it  merely 
as  a  change  in  velocity. 


For  the  pulse  shape  is  of  major  importance  and,  for  a  rectangular  pulse,  the 

peak  forces  generated  in  the  body  are  a  function  only  of  the  pulse  magnitude  (  Ye  ),  ir¬ 
respective  of  pulse  length. 


(1)  Linear  model  with  damping. 

So  far  the  discussion  has  been  confined  to  modes  in  which  the  human  body  behaves  as  a 
"single  degree  of  freedom"  dynamic  system,  of  the  type  showii  in  figure  19.  In  practice, 
this  happens  to  be  a  good  approximation  of  the  human  body's  response. 
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In  the  case  of  a  sitting  subject  experiencing  a  positive  spinal  (eyeballs  down)  acceleration, 
for  example,  we  might  give  physical  reality  to  figure  19  by  the  following  identifications. 

The  mass  corresponds  to  the  mass  of  the  head  and 
upper  torso. 

The  spring  corresponds  to  the  subject's  spine. 

The  damper  corresponds  to  the  distributed  damping 
in  the  spine  and  associated  tissues. 

In  Reference  10  Stech  showed  that,  if  the  effective  spring  rate  (  )  was  obtained  from 

Cadaver  data,  this  postulate  could  be  used  to  calculate  the  frequency  of  a  dynamic  model 
representing  the  human  body.  The  result  agreed  with  other  measurements  made  with  live 
human  subjects.  However,  it  is  important  to  realize  that  such  identifications  cannot  always 
be  made  (the  transverse  case  for  example)  and  that,  strictly  speaking,  they  are  not  necessary 
to  our  purpose. 


Figure  19.  Single  degree  of  freedom,  linear  dynamic 
system. 


We  shall  now  consider  the  mathematics  of  a  single  degree  of  freedom  system  in  order  to 
determine  the  lumped  parameter  coefficients  for  particular  tolerance  curves  of  the  type 
shown  in  figui'e  18,  when  the  dynamic  model  is  linear.  The  assumption  of  linearity  is 
made  for  two  reasons.  First,  because  the  resulting  equations  are  very  much  easier  to  use 
use  and  secondly,  because  available  experimental  data  is  not  precise  enough  to  permit 
identification  of  nonlinear  effects. 

We  define  the  following  linear  coefficients: 

Frequency  ,  / 
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Damping 


The  peak  force  in  the  spring  is  Dividing  by  the  mass  we  have  which  is 

the  "Dynamic  Response  Index.  " 


Note  that 


- 


peak  force 
effective  mass. 

We  assume  that  the  size  of  the  structure  (such  as  the  spine)  in  which  the  force  is  reacted 
is  proportional  to  the  mass  of  the  system  and  therefore  to  the  weight  of  the  subject.  Thus, 


^  proportional  to  peak  force 

structural  area 

/.  is  proportional  to  peak  stress. 

We  therefore  assume  that  the  Dynamic  Response  Index  (DRI)  is  a  measure  of  the  maximum 
stress  generated  in  the  body.  Thus,  the  DRI  is  the  most  logical  parameter  against  which 
to  correlate  injury  and  other  stress-induced  physiological  effects. 

For  zero  initial  conditions  ^^<J]i^is  given  by  the  following  equations: 

Impulsive  velocity  change 


(170) 


Short  period  acceleration  Yc 


(171) 


Now  we  have  assumed  that  any  line  in  a  tolerance  graph  is  defined  by  constant. 

Thus,  we  can^equate  equation  170  and  171  for  any  particular  tolerance  line,  with  values 
and  C  respectively. 
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The  "corner  duration"  in  figure  18  can  also  be  obtained  from  equation  172,  since 


(173) 


at  this  point.  Substituting  for  in  equation  172 


/  -h  e 


/If, 


This  equation  is  plotted  in  figure  20. 


(174) 


We  see  that  for  a  conventional  tolerance  graph,  one  parameter,  the  "corner  duration" 
defines  the  frequency  of  the  model  if  the  damping  ratio  is  known.  However,  it  is 
impossible  to  deduce  the  damping  ratio  from  tolerance  data.  This  must  be  obtained 
from  other  classes  of  experiment. 


b.  The  use  of  impedance  data  to  define  damping  and  nonlinearity. 

Probably  the  simplest  —  and  the  most  accurate  —  experiments  which  can  be  performed 
on  the  human  bod3' ;;  re  thos^  in  which  sinusoidal  vibratiou  is  imposed  upon  a  live  subject 
and  measurements  are  made  of  the  "Mechanical  Impedance"  defined  as 

Mechanical  Impedance  Jzf  -  - — 

Velocity  of  input  point 


A  second  common  experimental  method  involves  determining  Amplitude  Transmissability 
by  measuring  the  amplitude  of  oscillation  of  various  parts  of  the  test  subject’s  body  and 
dividing  these  readings  by  the  amplitude  of  the  vibration  input.  This  is  much  less  accurate 
than  impedance  measurements  because  we  have  to  assume  that  transducers  mounted  on  a 
subject  give  an  accurate  transcription  of  the  motion  of  that  part  of  the  body  to  v/hich  they 
are  attached.  In  practice  this  is  usually  impossible  for  frequencies  in  excess  of  50  c.  p.  s. , 
as  shown  by  the  work  of  von  Gierke  (Reference  4,  p.  155),  Significant  errors  may  occur 
at  much  lower  frequencies  if  the  transducer  attachments  are  not  carefully  designed.  Never¬ 
theless,  amplitude  ratio  measurements  are  qualitatively  useful  because  they  yield  a  much 
clearer  physical  picture  of  the  body's  modes  of  oscillation  and  thus  aid  in  a  clearer  under¬ 
standing  of  the  mechanism  involved. 


Figure  20.  Variation  of  nondimensionnl  comer  duration 

with  damping  ratio  for  a  recta ngiilur  acceleration 
pulse. 

ns 


As  variants  of  amplitude  transmissability,  some  experimenters  have  used  acceleration 
ratio.  Velocity  ratio  is  also  a  possibility.  Both  of  these  ratios  are  proportional  to  am¬ 
plitude  ratio,  so  long  as  the  motion  is  sinusoidal.  In  nonlinear  systems  however,  there 
may  be  marked  differences  between  them ,  so  that  it  is  important  to  define  the  actual  ratio 
being  employed. 


Figure  21.  Mechanical  impedance  of  a  linear,  single  degree 
of  freedom  dynamic  system  excited  at  the  spring 
base. 


The  mechanical  impedance  of  a  linear  dynamic  system  is  sketched  in  figure  21,  and  the 
peak  value  of  /  Z/  occurs  at  a  frequency  >  given  by 


"7^  ” 


(175) 


and  the  peak  value  of/2/  is  approximately 


*  -f- 


(176) 


These  relationships  are  derived  in  Ajjpendix  II. 


Since  it  is  not  necessary  to  instrument  the  test  subject  in  order  to  obtain  mechanical  im¬ 
pedance  readings,  this  test  procedure  is  the  most  precise  and  at  the  same  time  the  most 
satisfying  from  a  scientific  point  of  view.  It  should  be  noted,  however,  that  resonant 
frequencies  caneot  be  obtained  directly  from  impedance  data.  If  there  is  appreeiabU' 
damping  in  the  syotern,  the  frequencies  at  which  maximum  impedance  occurs  are  con¬ 
siderably  higher  than  the  resonant  frequency.  The  relationship  between  them  is  (riven 
in  figure  22.  The  maximum  impedance  for  a  linear  system  is  plotted  in  figures  2;i  and 
24. 

An  impedance  curve  can  tell  us  the  damping  in  a  system.  From  figure  21  it  is  obvious 
that  we  know / Z /  and  ,  the  frequency  at  which) Z/ ^4* occurs.  Thus,  from 

figure  24  we  can  determine  the  damping  ratio  ?  ,  if  we  know  the  effective  mass. 

It  should  be  noted  that  we  are  almost  always  concerned  with  the  parameter  /z^ ,  not 
lZ\  alone.  It  is  strongly  recommended  therefore  that  future  impedance  measurements 
be  reported  in  the  form  of  ,  rather  than  (X.(  .  The  former  is  a  more  fundamental 

parameter  and  should  reduce  scatter  by  eliminating  the  effect  of  weight  variation  between 
subjects.  Moreover,  the  variability  of  /2//jk  can  be  statistically  analyzed,  whereas  a 
knowledge  of  the  /zy  distribution  is  meaningless  because  it  may  reflect  nothing  more 
than  a  variation  in  the  weights  of  the  subject  tested. 

Once  the  damping  factor  is  known,  the  undamped  natural  frequency  can  be  obtained  from 
figure  22. 

If  amplitude  transmission  measurements  are  made,  the  relationship  betw'een  the  resonant 
frequency  and  the  peak  impedance  frequency  will  also  enable  the  damping 

ratio  to  be  read  from  figure  22. 

Impedance  measurements  with  the  human  body  reported  by  Coermann  in  Reference  f) 
enable  preliminary  estimates  of  the  damping  coefficient  "c  to  be  made  for  relatively 
small  amplitude  oscillations,  at  frequencies  up  to  20  cps.  The  values  derived  by 
Coermann  are  not  as  accurate  as  they  could  bo,  however,  because  he  has  used  an  approxi¬ 
mate  equation  for  maximum  impedance.  In  addition,  he  used  only  the  peak  impedance  to 
determine  c"  ,  whereas  the  ratio  of  the  resonant  frequency  to  the  maximum  impedance 
frequency  can  also  be  used  to  obtain  a  second  reading. 

For  one  subject  (code  R.  C.)  Coermann’s  data  can  be  summarized  as  follows,  taking  the 
subject  mass  as  185/32.2  -  5.  75  slugs. 


Sitting  Erect  ^ 

(Coermann)  (Latham) 

Sitting 

Relaxed 

Standing 

Erect 

Resonant  frequency 
(cp-  ^ 

5.  2 

5.5 

4.65 

5.05 

Frequency  for 

6.  3 

5.  3 

5.  9 

^  calculated 
by  Coermann 

0.285 

0.  325 

0.  37 

/  mAk 

(d3Tie  X  sec/cm  x  10®) 

6.65 

5.08 

5.  2 

Z/  (lb  sec /ft) 

465 

348 

356 

2.005 

1.82 

1.672 

based  on  /z/^^jr 

0.281 

0.317 

0.  352 

1.21 

1.14 

1.168 

based  on 

0.352 

0.305 

0.  325 

Mean  value  for 

0.306 

0.316 

0.  349 

Natural  frequency 

6.1 

5.14 

5.67 

t^/ZfT  (cps) 


For  practical  purposes  we  may  neglect  the  variation  betw<,cn  sitting  erect  and  relaxed. 
Thus,  the  final  figures  become 

Damping  factor  in  sitting  position,  ^  =  0.31 
Damping  factor  standing  stiff-legged,  C  =  o.  35 

c.  The  application  of  rebound  data  to  determine  dynamic  constants. 

Impedance  measurements  enable  us  to  estimate  the  damping  in  the  humrn  body  if  we  are 
able  to  excite  the  mode  which  is  of  importance.  We  cannot  obtain  the  full  deflection  of  the 
body’s  structure  that  would  be  experienced  in  an  impact  test,  for  example.  It  is  obvious 
that  the  mean  effective  damping  in  an  impact  might  be  different  from  that  nieasured  in  a 
relatively  low  amplitude  impedance  test. 
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-fl.#  =  FREQUENCY  FOR  MAX.  IMPEDANCE 
=  UNDAMPED  NATURAL  FREQUENCY 


RESONANT  FREQUENCY 


DAMPING  COEFFICIENT  RATIO  »  C/t0i 


Figure  22.  Variation  of  maximum  impedance  frequencv 
ratio  with  damping  coefficient,  for  a  linear 
system. 
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i^’lgure  23,  variation  of  maximum  impedance  with  (lamping 

ratio,  in  terms  of  the  undamped  natural  frequency  oH 
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Figure  24,  Variation  of  maximum  inpedance  with  damping 

ratio,  in  terms  of  frequency  for  peak  impedance. 
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Secondly,  impedance  measurements  can  bo  used  to  determine  the  degree  of  nonlinearity 
only  if  a  suitable  centrifuge  is  available  for  use  in  the  experiments. 

Rebound  testing  appears  to  offer  a  means  of  circumventing  both  these  limitations. 

(1)  Rebound  of  a  linear  system. 

If  a  linear  dynamic  system  is  subjected  to  an  impulsive  velocity  change,  such  as  l)y  the 
drop  test  indicated  in  figure  25 ,  it  will  rebound  away  from  the  free  surface  on  which  it 
impacted. 


Figure  25*  Impact  case  for  a  simple  dynamic  system. 


The  equation  of  motion  is.  of  course, 

#  •  •  X  ^ 


the  initial  condition  being  and  the  solution 

(178) 

(179) 

Where  f  =  ^  ^  c/c^ 

"Take  off"  occurs  either  when  the  total  spring  and  damper  force  falls  to  zero,  or  when 
the  spring  deflection  is  zero.  For  the  first  of  these 


«  o 


(180) 


at  "take  off". 
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Substituting  equation  178  and  179  for  S  and  ^  equation  180  gives  the  result 

XTr” —  ^ Zc^  s  ^(say) 


(181) 


Substituting 


into  equation  179,  the  take  off  velocity  is 

^  5-  l^-r-9  V  <fl 


(182) 


Equation  182  is  plotted  in  figure  28,  together  with  the  velocity  which  exists  when  the 
spring  deflection  is  zero. 


For  the  case 


S  =0, 


=  0,  so  that 


\t  *  ^TT  ^  ir  * 

~  Z UMA  If 

^  -e  ^ 


Mrr 

A 

4 


Figure  26.  Variation  of  rebound  velocity  with 
damping  in  a  linear  system. 


(183) 
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Consequently,  if  we  measure  the  rebound  velocity  of  live  human  subjects,  we  can  immedi¬ 
ately  determine  the  equivalent  linear  damping  from  figure  26  or  equation  183,  without  any 
reference  to  the  frequency  of  the  system. 


Once  we  know  the  damping,  the  time  to  take  off  (  ^  )  gives  us  the  frequency  of  the  systeni, 
since  from  equation  181 


^to 


(184) 


This  relationship  is  plotted  in  figure  27* 

Thus,  we  see  that  rebound  measurements  can  give  us  a  direct  and  precise  indication  of 
both  damping  and  frequency,  under  the  type  of  loading  history  with  which  we  are  most 
concerned  in  body  dynamics.  '  "he  theory  given  above  is  for  true  "bounce  testing"  where 
the  subject  is  dropped  onto  a  rigid  surface.  In  Reference  12  Hirsch  reports  tests  in  which 
take  off  is  caused  by  vertically  accelerating,  and  then  decelerating,  a  platform  on  which 
the  subject  is  standing.  This  data  is  more  difficult  to  analyze  because  the  acceleration 
time  history  occupies  a  finite  time.  Some  of  Hirsch's  data  is  plotted  in  figure  28,  as  are 
also  boundaries  for  an  undamped  single  degree  of  freedom  systeni  driven  by  rectangular 
pulses.  It  seems  obvious  from  the  relatively  small  scatter  in  Hirsch's  data  that  good 
estimates  of  both  damping  and  frequency  could  be  made  by  using  an  analog  to  reduce 
this  data.  Also,  if  such  tests  were  carried  out  with  a  carefully  regularized  acceleration¬ 
time  profile,  the  damping  and  frequency  of  the  subjects  could  be  obtained  by  a  simple 
analysis  of  the  type  developed  earlier  in  this  section. 

(2)  Rebound  of  a  nonlinear  system. 

If  the  apparent  frequency  of  the  human  body  increases  with  impact  velocity,  during  re¬ 
bound  testing  of  the  type  described  in  the  previous  section,  we  may  feel  reasonably  con¬ 
fident  that  this  is  due  to  nonlinearities  in  the  system,  and  most  probably  to  a  progressive 
increase  of  spring  rate  with  deflection. 

It  should  be  noted  that  this  type  of  experimental  work  offers  the  hope  of  a  major  breakthrough 
in  the  area  of  body  dynamics.  Consequently  it  is  recommended  that  a  suitable  experimental 
program  be  implemented  with  all  possible  priority. 

d.  Some  examples  of  the  analysis  of  existing  data. 

At  the  time  of  writii^  this  report,  insufficient  data  are  available  to  make  a  detailed 
analysis  of  the  dynamic  models  which  best  describe  the  dynamic  behavior  of  the  human 
body.  The  process  of  collecting  the  available  data  is  still  far  from  complete.  However, 
it  is  of  interest  to  construct  some  provisional  models  in  order  to  illustrate  the  use  of  the 
techniques  described  in  the  earlier  sections. 


NONDIMENSIONAL  TIME,  AT 


Figure  27.  Effect  of  model  damping  on  time  to  7.ero  force 

i  \  and  v.firo  extension  f  ^  fnt*  an  immilHivo 
input  to  a  damped,  linear  sybtem. 
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Undamped  dynamic 
^  T  system  v;ith  positiv 
I  puLse  only. 


Reference  12  data  p  ants 
for  human  subjects,  with 
arbitrary  negative 
pulse. 


Undamped  dynamic 
system  with  infinite 
deceleration  after 
pulse. 


DURATION  OF  POSITIVE  ACCELERATION  PULSE 


UNDAMPED  NATURAL  PERIOD  OF  DYNAMIC  MODEL 


Figure  28.  Take  off  velocity  as  a  function  of  pulse  duration 
for  an  undamped  dynamic  system ,  compared  with 
Hirsch’s^^  data  for  a  standing  man. 
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(1)  Spinal  model. 


FojT  tliC  GasG  of  positiVG  spinal  accsloi^ation,  it  lias  l#GGn  fonncl  ttiat  IKg  availalilB  data  on 
spinal  injury  support  the  hypothesis  that  a  single  degree  of  freedom  nonlinear  damped 
system  is  a  good  representation  of  t’  e  human  body.  An  estimate  of  the  nonlinearity  of 
the  system  has  been  made  in  Referet  ce  10. 

The  mean  age  for  U.  S.  Air  Force  aircrew  personnel  is  27.  9  years.  This  age  is  used 
as  a  basis  for  estimating  tolerances  to  acceleration. 

Using  data  from  experiments  with  cadaver  vertebrae,  and  a  knowledge  of  the  variation 
of  their  strength  characteristics  with  age  (Reference  11),  it  can  be  shown  that  the  steady- 
state  acceleration  level  corresponding  to  50%  probability  of  injury  for  a  subject  aged  27.  9 
years  is  approximately  21.  3g. 

From  Table  1  of  Reference  11,  corresponding  to  these  conditions  we  have 

=  0. 2245 
=  8. 42  cps 

'Ait 

Figure  20  then  shows  that 

=  2.27 

Since  ,  the  real  time  corner  duration  is  therefore 

Atf  ~  — •  r^~'~ —  =  0.0428  secs. 


.  ^ 


(185) 


8.42  X  2 


(186) 


The  d5mamic  overshoot  of  a  system  subjected  to  a  zero  rise  time  acceleration  input 
and  having  a  damping  coefficient  of  ^  =  0.2245  is,  from  equation  171, 


•  « 

>'c 


=  1.485 


(187) 


This  means  that  the  critical  impulsive  velocity  change  for  50%  probability  of  injury  is 


_  21,3  X  32.2 


1.485 


X  0.0428 


=  19.  75  ft/sec 


(188) 


Using  the  values  expressed  in  equations  186  through  188,  we  are  able  to  define  the  50% 
probability  of  injury  curve  for  the  case  of  a  rectangular  pulse  acceleration  input.  This 
is  shown  in  figure 2 9. 


65 


ACCELERATION  MAGNITUDE 


NOTE:  This  injury  curve  is  based  upon  dynamic  models  of  the 

human  body  using  available  e3q>erimcntal  data.  It  is 
subject  to  revision  as  additional  data  becomes  available. 

The  only  tolerance  curves  authorized  by  government 
procurement  agencies  for  the  use  as  design  limits  are 
those  defined  in  HIAD  or  other  applicable  procurement 
specifications. 


0.01  0.0428  0.1  1.0 

DURATION  OF  ACCELERATION  PULSE  N  SECS. 


Figure  29.  50%  probability  of  injury  curve  for  a  r^’ctangula 

acceleration  pulse  in  li.e  spinal  direction. 


ac 


Summarizing  the  model,  we  have 


=  52. 9  rad /sec 

“c 

=  0.2245 

Corner  duration 

=  0.428  see 

Critical  velocity  cliange  ~  19.  75  ft/sec  (A't  K. 

Critical  zero  rise  time  v 

acceleration  =  14.55  g 

It  should  be  noted,  as  pointed  out  in  Reference  11,  that  head  involvement  occurs  when 
the  duration  of  the  acceleration  pulse  becomes  very  short.  At  the  present  time,  insuf¬ 
ficient  data  are  available  to  determine  the  limitations  associated  with  this  degree  of 
freedom,  however. 


(2)  Transverse  model  (soft  head  restraint). 


There  is  a  substantial  amount  of  data  to  indicate  that,  for  the  positive  transverse  case 
the  50%  injury  level  is  given  by 

=  53. 6  ft/sec 

'i*,r  ■ 


for  rectangular  acceleration  pulses  with  zero  rise  time,  and  that  a  single  degree  of 
freedom  dynamic  model  represents  an  adequate  description.  These  figures  will  probably 
change  somewhat  as  the  process  of  data  analysis  continues,  however. 


From  equation  173 

Thus,  if  we  know  the  damping  ratio  C  ,  we  could  obtain  the  frequency  from  figure  20. 
For  the  probable  range  of  values  we  have 


f*Atr 


53,0 


40.0  X  32.2 


=  .0416  secs. 


0 

0.1 

0.2 

0.3 

0.4 

*» 

2.0 

2.6 

3.45 

4.72 

6.83 

<0  = 

48. 1 

62.5 

82.9 

113.4 

164. 1  rads /sec 

Until  rebound  tests  have  been  carried  out,  the  only  source  of  information  on  damping  in  the 
transverse  mode  is  due  to  a  single  impedance  curve  presented  by  Coermann  in  Reference  8^ 
This  gives  a  peak  impedance  of  432  lb  sec /ft  at  a  frequency  -  7.  6  cps.  Assuming 

tiiat  the  weight  of  the  subject  was  152  lbs.  , 


1.  92 


and  from  figure  24. 


^  =  0.298 


Since  this  result  is  for  a  "semisupine"  position,  and  is  a  single  data  point,  not  much  weight 
can  be  attached  to  this  result.  It  is  extremely  interesting  to  note  that  the  peak  impedance 

frequency  is  47.  7  rads /sec.  ,  however,  since  this  is  comparable  with  the  values  deduced 
from  pulse  tolerance  data. 
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3.  DYNAMIC  CONSIDERATIONS  IN  THE  DESIGN  OF  EXPERIMENTS 
WITH  LIVE  HUMAN  SUBJECTS. 

Dynamic  models  are  of  considerable  value  in  the  design  of  physiological  experiments, 
and  can  save  a  great  deal  of  time  and  money.  For  very  short  acceleration  durations , 
where  velocity  change  only  is  of  importance,  there  is  little  point  in  using  a  sled  or 
HYGE  accelerator,  since  exactly  the  same  results  an  be  achieved  by  drop  testing. 

The  size  of  test  equipment  is  also  an  important  factor.  Cases  have  arisen  in  the  past 
where  short  period  inputs  have  been  attempted  on  equipment  whose  maximum  working  stroke 
automatically  limits  it  to  "impact"  or  impulsive  velocity  change  experiments.  Foi  .unately, 
the  dynamic  models  so  far  established  allow  us  to  define  test  equipment  parameters  quite 
precisely.  The  theoretical  considerations  involved  are  examined  in  the  next  section. 

a.  Test  equipment  requirements  for  determining  the  tolerance  limit  of  a  single 
degree  of  freedom  dynamic  model. 


Figure  30.  Rectangular  pulse  tolerar  ce  graph  for  a  single 
degree  of  freedom  linear  dynamic  system  with 
zero  damping. 


Consider  the  nondimens ional  tolerance  graph  sketched  in  figure  30  for  an  undamped  model. 
The  velocity  change  associated  with  the  rectangular  pulse  input  upon  which  this  graph  is 
based  will  be 


.... 


(189) 
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A  suitable  nondimensionalization  for  is  therefore 


•  • 


(190) 


Equation  190  is  sketched  in  figure  31. 


1.0  a.  0  3.  0 

NONDIMENSIONAI.  DURATION  AT  «•  i^At 


Figure  31.  Nondimensional  velocity  change  required  in 
testing  with  a  rectangular  pulse. 

It  should  be  noted  that  true  impace  velocity  changes  occur  for  durations  oi  AT ^  1.0, 
and  true  short  period  accelerations  in  the  range  of  AT  >  3.  0.  A  velocity  change  of  at 
least 

is  required  in  the  short  period  range. 
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The  associated  stroke  required  over  the  acceleration  period  is  given  by 


(191) 


so  that  a  suitable  nondimensionalization  is 


(192) 


This  relationship  is  shown  in  figure  32. 


0  0.1  0.2  0.  3  0.4 


NONDIMENSIONAL  DURATION 

Figure  32.  Variation  of  nondimensional  accelerator  stroke 
with  acceleration  duration. 


Thus  we  can  define  the  acceleration  stroke  required  as  follows; 

S’  <  { 


iO 


s,  >  I 


HAff 


Impact 

Short-period 
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(193) 


We  have  seen  that  the  parameter  which  most  conveniently  discriminates  between  impact 
and  short  period  acceleration  is  iJT  ,  the  nondimen sional  acceleration  duration. 


For  impact  accelerations  K.  I 

For  short -period  accelerations  4*2^  >  3 


(194) 


The  stroke  of  an  accelerator  which  produces  a  rectangular  pulse  is 


s,  ■=  = 

(195) 

- 

(196) 

X  ^ 

For  impact  experiments,  ^ 

2 

•  * 

(197) 

For  short  period  experiments  5^  ^ 

X  ^ 

(198) 

Equations  197  and  198  agree  with  equation  193,  as  we  should  expect.  The  actual  test  rig 
stroke  limitations  involved  are  plotted  in  figures  33  and  34  for  spinal  acceleration,  and 
figures  35  and  36  for  transverse,  using  the  existing  dynamic  models  (Reference  2).  The 
operating  range  of  S,  in  figures  33  through  36  is  the  range  between  the  vertical  cross- 
latched  lines.  Limited  sample  testing  at  levels  of  less  than  0. 1%  probability  of  injury  is 
not  considered  worthwhile.  On  the  other  hand,  testing  at  levels  oi  greater  than  10%  proba¬ 
bility  of  injury  is  undesirable,  from  safety  considerations.  We  may  summarize  these  re¬ 
sults  as  in  Table  1  below: 


TABLE  1 

Test  Rig  Stroke  Limitations 

Impact  (impulsive  velocity  change)  Testing 

Stroke  must  be  less  than  ....  1. 29  in.  for  spinal  (0. 1%  injury) 

2. 24  in.  for  spinal  (50%  injury) 

3. 45  in.  for  transverse  (0. 1%  injury) 
6. 70  in.  for  transverse  (50%  injury) 


Short-period  A  cceie ration  Testing 


Stroke  must  be  greater  than  .  3.  80  in. 

6.61  in. 
8. 63  in. 
15.  02  in. 


for  spinal  (0. 1%  injury) 
for  spinal  (50%  injury) 
for  transverse  (0. 1%  injury) 
for  transverse  (50%  injury) 
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PROBABILITY  OF  INJURY  MAXIMUM  ACCELERATION  STROKE  (INCHES) 


0.10 


0.01 

0.001 


0.  01 


0. 10 


Impact  test  rig  strokes  for  spinrl  cxccier:  - 
tion  of  loss  than  0.04  sec,  diir^’thm,  u  ing 
live  humrn  subjei'ts. 


Figure  33, 


PROBABILITY  OF  INJURY  MINIMUM  ACCELERATION  STROKE  (INCHES) 


10.0 


1.0 


0.10 
0.  001 


0.01 


0.10 


Figure  34.  Short  period  test  rig  stroke  for  spinal 
acceierution  of  gi’eater  than.  0,  00  sec, 
duration,  using  live  human  subjects. 
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PROBABILITY  OF  INJURY  MAXIMUM  ACCELERATION  STROKE  (INCHES) 


PROBABILITY  OF  INJURY  MINIMUl^I  ACCELERATION  STROKE  (INCHES) 


Figure  36-  Short  period  test  rig  stroke  for  transverse 
testing  with  live  human  subjects. 
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These  results  are  particularly  significant,  since  they  indicate  much  more  stringent  test 
rig  limitations  than  generally  realized.  In  the  spinal  case,  for  example,  the  structure 
supporting  the  test  subject  should  be  extremely  rigid,  and  absolutely  no  cushioning  ma¬ 
terial  can  be  used  if  the  results  are  to  accurately  portray  the  response  of  the  human  body 
alone.  *  Also  of  interest  is  the  observation  that  a  very  limited  stroke  device  —  such  as 
a  small  diameter  HYGE  tester  —  cannot  be  used  for  short-period  experiments.  Equally 
evident  is  the  fact  that  "dirt  drops,  "  such  as  those  carried  out  by  Holcomb  with  the  B-58 
capsule,  and  other  tests  involving  some  resilience  in  either  the  test  rig  or  the  surface  upon 
which  it  imnacts,  cannot  be  reduced  to  give  useful  data,  unless  the  acceleration-time  history 
is  accurate  y  recorded.  Usually  it  is  not. 

In  the  case  of  transverse  experiments,  the  reverse  is  true.  Almost  any  practical  support 
and  restraint  system  may  be  used  for  impact  testing,  so  long  as  the  deflection  under  40  g 
static  load  does  not  exceed  about  two  inches.  Correspondingly,  larger  strokes  are  required 
for  short-period  investigations  and  these  can  probably  only  be  carried  out  on  specialized 
rocket  sled  and  similar  facilities. 


Thus  the  type  of  test  rigs  needed  may  be  roughly  summarized  as  follows: 


Impact 


Short- Period 


Spinal  Drop-Test  (very 

stiff  support) 


HYGE  Tester  (very 
stiff  support) 


Transverse  Drop-Test 


Rocket  Sled,  Daisy 
Track,  Ejection 
Tower,  Etc. 


For  practical  purposes,  the  requirements  in  lateral  testing  will  approximate  those  for 
transverse. 


b.  The  influence  of  the  restraint  system. 

A  restraint  system  is  some  form  of  resiliency  between  a  human  body  and  the  structure  of 
the  vehicle  or  vehicle  subassembly  with  which  it  is  associated.  When  the  vehicle  is  sub¬ 
jected  to  an  acceleration,  the  existence  of  the  restraint  system  means  that  the  acceleration¬ 
time  history  felt  by  the  man  is  not  the  same  as  for  the  vehicle.  This  modification  of  the 
vehicle's  acceleration  can  be  favorable  or  unfavorable,  depending  upon  the  dynamic  char¬ 
acteristics  of  the  restraint  system. 


♦  It  is  reported  that  Swearingen^  used  a  hydraulic  snubber  in  many  of  his  "impact" 
experiments,  so  that  they  cannot  necessarily  be  regarded  as  "impact"  experiments. 
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A  good  restraint  system  will  reduce  the  physiological  effect  of  an  acceleration  because 


(a)  It  increases  the  effective  stroke  through  which  the 
man  is  accelerated  to  a  given  velocity  change,  and 

(b)  It  increases  the  rise  time  of  the  acceleration,  thus 
reducing  the  dynamic  overshoot  of  the  man. 

A  bad  restraint  system  will  increase  the  physiological  effect  of  an  acceleration  because 

(a)  It  permits  the  man  and  the  vehicle  to  develop  relative 
motion  before  it  bottoms  out,  thus  increasing  the  shock 
effect  o'  a  given  acceleration  input. 

The  detailed  effects  of  restraint  systems ,  and  the  way  in  which  they  modify  the  physiological 
effect  of  inpu  accelerations  have  been  discussed  in  Part  I  of  this  report  and  can,  in  principle, 
at  least,  be  c  nsidered  in  the  analysis  of  test  data.  However,  their  presence  is  a  complicat¬ 
ing  factor.  At  the  present  time,  it  seems  logical  to  say  that  absolutely  no  series  resiliency 
should  be  permitted  in  spinal  tests,  and  that  the  maximum  restraint  deflection  in  lateral  and 
transverse  tests  should  be  limited  as  much  as  possible,  preferably  to  less  than  one  inch,  and 
certainly  to  less  than  two  inches. 

When  the  use  of  a  restraint  resiliency  is  unavoidable,  its  dynamic  characteristics  should  be 
recorded  in  as  much  detail  as  possible,  so  that  its  effects  can  be  calculated  and  the  basic 
input  pulse  shape  to  the  human  body  determined  later  by  a  device  such  as  the  Frost  Restraint 
Analyzer. 
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APPENDIX  I 


SOLUTIONS  OF  THE  SINGLE  DEGREE  OF  FREEDOM  LINEAR 
DIFFERENTIAL  EQUATION  WITH  CONSTANT  ACCELERATION  INPUT 


The  basic  equation  in  linear  dynamics  is 

j  +  ui  j’?  -i/c 

#  • 

where  is  some  function  of  time. 

1.  The  Complementary  Function 

When  the  damping  is  subcritical,  that  is  to  say,  when  C  <  ^  the  equation 

$  -f  Qci  i-  0^^  S 


1(1) 


1(2) 


has  the  solution 


u> 


1(3) 


or 


,  if  r  =  9^5 


S  =  J  I  coYSK  ^  -i-WSll  } 


1(4) 


where  (Si  and  (Si  are  initial  conditions  and 

This  solution  is  the  Complementary  Function  in  the  solution  of  Equation  1(1). 

2.  The  Particular  Integral 

The  solution  to  Equation  1(1)  contains  a  Particular  Integral  determined  by  the  function 

/(f).  ^ 

This  is  given  by  ^ ^  2)  +  ^  ^ 


where  m  denotes  the  operator 
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3.  Constant  Acceleration  Input  ^ _ With  Sub-Critical  Damping. 

•  • 

In  the  special  case  when  it)  -  K  =  a  constant, 


y 

tc* 


because  7>  (  )  and  J!>  (  u  )  =  0,  must  be  added  to  the  complementary 

function  S  oi  Equation  1(4)  noting  that  the  initial  displacement  in  Equation  1(4)  is  now 
instead  of  (  S  L  •  Equation  1(4)  now  becomes 


where 


f,  ’ 

z  = 


1(6) 


1(7) 


"^•/a/i  -tZf,+  f,‘ 

f Vr- 

<^/u) 


7/  ~  2  *  * 


£0 


If  we  write 


fK  ^  t  -  tiJP* 

Yc 


^  T  ^  i 


1(8) 


The  displacement  S  has  maximum  and  minimum  values  when  the  displacement  velocity 
J  =  0.  These  are  given  when 


€ 


5>)  ^(A  \  ^cuiS)  «  o 
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That  is,  when 


1(9) 


where 


Equation  1(9)  is  satisfied  when 


X,  =  •<ir  where  ^  =0,  1,  2,  3.... 

mA> 


1(10) 


Now  from  Equation  1(8), 


—  '  "  f  %  • 


1(11) 


1(12) 
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=  C-if  fjn. A 

jnfvT7'  /-2<r*f 

f  ^ 


/^f 

Substituting  in  Equation  1(7), 


Jf-oV/f-  ^  ft’- \ 


1(13) 


>i 


where  ACr  ' 


•/ 


^t-  2^^  -f 


1(14) 


1(15) 


Note  that,  for  zero  initial  values  of  S  and  S  Equation  1(7)  becomes 

-?*j  t 

7. 


1(16) 


where 


S«t*  ^  ^  —  "c 


The  maximum  value  of  is  then 

X 

-f2? 


1(17) 


and  occurs  when  ^ 


w  -rr- 

"A 
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APPENDIX  n 


THEORY  OF  MECHANICAL  IMPEDANCE  FOR  A  SINGLE 


DEGREE  OF  FREEDOM,  LINEAR  DYNAMIC  SYSTEM 


In  this  Appendix, we  consider  a  damped,  linear  system,  for  which  the  equation  of  motion 
is 


11(1) 


Writing  equation  11(1)  in  Laplacian  form,  and  neglecting  transients, 


or  ; 


The  inverse  of  +his  is 


n(2) 


i. 


M 


Neglecting  transients,  this  reduces  to 


n(3) 


where 


n(4) 


Now  if  A  sin  at  is  the  input  acceleration,  the  input  velocity  is 

A  A^Jtt 


n(5) 
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The  force  in  the  system  is 

-i-  h 

From  equation  11(3) 

*  JL  ee^  ^,) 

^  -  r-  - - 


From  equations  n(3)  and  n(5)  the  force  in  the  system  is  given  by 

P  ^  Am.  +'^>y  ■*  (-o-t  f 7t 


where 


and 


f,- 


Now  from  11(5) ,  we  have  that 


Jc  ^ 


n(6) 


11(7) 


n(8) 

n(9) 

n(io) 


n(ii) 


Dividing  equation  11(7)  by  11(11)  the  mechanical  impedance  is 

5^  (At  t  y,  -t^,) 

Writing  ^  ^  and  ^  -  2c  ,  the  modulus  is 


n(12) 
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n(i3) 


•B  ^ 


Obviously  the  phase  angle  is 


n(14) 


Of  interest  is  the  value  of  / Zj  when  ,  so  that  =  1 


n(i5) 


In  general  the  equation  for  is  of  more  value  than  the  equation  for  lz/u»  , 

since  there  is  usually  no  way  of  knowing  the  correct  value  for  to 


From  equation  n(  13), 


m  ~  h 


A’'* 


n(i6) 


Differentiating  with  respect  to  jo  , 


—  /'III.)  « 

[tutoj 


■k)CI’* 


where  ^  ^  —  JL 

Equating  to  zero  to  find  the  value  for  /  Z 


0 


4"^^  4  /if f 
(  -  VVc‘-  2) 


n(i7) 


11(18) 


n(i9) 
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Note  that  as  ?  0,  1. 0 

as  ?■  — »  0.77,  ^  -►ott 

Another  limit  of  Interest  is  /  Zj  .  „  as 


•Mtt 


0. 


From  equation  11(16),  as  j9 


1.0 


From  equation  n(16)  note  that  as  )- 

'Zc 

A  damped  system  oscillates  at  a  frequency 


X 


.-I  — 


11(20) 


n(2i) 


Thus  when  a  measurement  of  the  frequency  ^  is  available,  as  well  as  an  impedance 
measurement,  we  may  determine  ?  from  the  ratio 


<A>  ^  ^  n(22) 


In  many  cases  the  frequency  for  maximum  amplitude  of  the  test  subject  (resonance)  is 
also  determined. 


Since  S  =  —  ('ll 

I  =  y'-*  i 

The  transmission  factors  (output/input)  are  therefore 


n(23) 
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1 


Displacement : 


Velocity; 


Acceleration: 


/-  K 


Tj-  =  /  -  f/p. 

Tj.  =  /  - 


From  equation  11(3), 


Also 


s>‘ 

s-  /4  ^:dr5-A‘^ 

4  «  -  A  ^  il'i- 
XI 

*4  S  C  —  4.  eoiSlt 

/It. 


n(24) 


n(25) 


n(26) 


Substituting  equations  11(25)  and  a(26)  into  11(24),  the  modulus  of  the  transmission  factor 
becomes 


that  is,  the  transmission  factor  is  the  same  for  acceleration,  velocity  or  displacement, 
so  long  as  the  motion  is  sinusoidal. 

We  find  the  "resonant  frequency"  for  which  “T  is  a  maximum  by  differentiating  equation 
11(27),  equating  to  zero,  and  solving  for 

This  gives  us  ^  — '  - 

toj  y7-  2f* 

In  other  words ,  the  resonant  frequency  is  higher  than  the  undamped  natural  frequency  cO  , 
although  the  free  oscillation  frequency  ^  is  lower. 
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